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ENGINEERING AEROTHERMAL ANALYSIS FOR 
X-34 THERMAL PROTECTION DESIGN 


Kathryn E. Wurster*, Christopher J. Riley + , and E. Vincent Zoby ++ 
NASA Langley Research Center, Hampton, VA 


Design of the thermal protection system for any hypersonic flight vehicle requires determina- 
tion of both the peak temperatures over the surface and the heating-rate history along the flight 
profile. In this paper, the process used to generate the aerothermal environments required for the X- 
34 Testbed Technology Demonstrator thermal protection system design is described as it has evolved 
from a relatively simplistic approach based on engineering methods applied to critical areas to one 
of detailed analyses over the entire vehicle. A brief description of the trajectory development lead- 
ing to the selection of the thermal protection system design trajectory is included. Comparisons of 
engineering heating predictions with wind-tunnel test data and with results obtained using a Navier- 
Stokes flowfield code and an inviscid/boundary layer method are shown. Good agreement is dem- 
onstrated among all these methods for both the ground-test condition and the peak heating flight 
condition. Finally, the detailed analysis using engineering methods to interpolate the surface-heat- 
ing-rate results from the inviscid/boundary layer method to predict the required thermal environ- 
ments is described and results presented. 


Introduction 

Continued human presence in space will soon re- 
quire upgrade or replacement of the current Space Trans- 
portation System. Joint NASA/industry partnerships have 
been formed to develop and test the technology required 
to provide cost-effective, yet reliable access to space. This 
is the primary goal of the Reusable Launch Vehicle 
(RLV) Technology Program. 1 Development and flight 
of the suborbital X-34 vehicle is an integral part of the 
RLV program and is intended to address numerous op- 
erations, manufacturing, and other technology issues that 
can best be evaluated through demonstration on an ac- 
tual flight vehicle. One of the key technologies to be dem- 
onstrated is that of an operationally effective thermal 
protection system (TPS). This paper describes the meth- 
odology by which the aerothermal environments for the 
X-34 are predicted in sufficient detail to allow design of 
the TPS such that the survivability, as well as the reus- 
ability, of the flight vehicle is ensured to a high degree 
of certainty. The primary focus of this paper is to de- 
scribe how wind-tunnel and computational -fluids results 
have been coupled with engineering methods to predict 
the time-dependent thermal environments required for 
the design of the X-34 TPS. 
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Aerospace vehicles, such as the X-34, that reach hy- 
personic speeds generally require protection of the ve- 
hicle structure from the aerodynamic heating at these con- 
ditions. TPS can be either of an active or a passive vari- 
ety. Active systems typically use onboard coolant pumped 
through “pipes” behind the external surface or employ 
heat pipes which use capillary action to redistribute the 
heat from hot areas of the vehicle to cooler areas. A pas- 
sive system, such as that employed on the X-34, takes 
advantage of the insulative properties of the TPS materi- 
als to hold a significant portion of the incident heat until 
it is radiated or convected away from the vehicle. The 
two critical parameters that determine the TPS options 
available for a given vehicle and trajectory are the peak- 
heating rate and the integrated heating over the time of 
the flight profile. The former determines the maximum 
temperature environment, and thus the material options, 
and the latter determines the thickness distribution of the 
insulative TPS material required to protect the surface. 
Except for the nosecap and wing leading edges (where 
Silicone Impregnated Reusable Ceramic Ablator 
(SIRCA) 2 tile is used), the X-34 is protected with flex- 
ible-blanket or felt systems. At the request of the lead 
contractor. Orbital Sciences Corporation (OSC), the Lan- 
gley Research Center (LaRC) effort concentrated on the 
surface areas to be protected by the flexible systems. The 
flexible blanket and felt systems are described in detail 
in Ref. 3. Preliminary estimates of the heating to the stag- 
nation areas will be included in this paper, but the bulk 
of the stagnation area analysis, both environment pre- 
diction and SIRCA tile design, has been performed at 
the Ames Research Center (ARC) and is described in 
Refs. 2 and 4. 



Coast 


In the past, engineering methods have been used al- 
most exclusively to predict external thermal environ- 
ments and define TPS requirements for hypersonic flight 
vehicles. The advent of high-speed computers that re- 
duce the time required to produce a computational solu- 
tion has enabled incorporation of detailed flowfield meth- 
ods much earlier in the design process. However, lim- 
ited resources and the complexity of the typical vehicle 
configuration generally preclude the generation of a suf- 
ficient number of computational solutions to adequately 
define the heating histories for TPS definition. Similarly, 
the recent development of experimental wind-tunnel 
techniques, 5 whereby heating distributions over the ve- 
hicle can be measured rapidly and at relatively low cost, 
also allows incorporation of experimental results at an 
early stage of the design process. These measured data 
are currently limited by the inability of the wind-tunnel 
to reproduce the actual flight environment (flight 
Reynolds numbers and reacting chemistry). Therefore, 
engineering techniques are used in this study to integrate 
the computational and the wind-tunnel results over the 
trajectory profile to provide the designers with sufficient 
information to adequately define the TPS requirements. 
The methodologies used to define the aerothermai envi- 
ronment and the process used in selection of trajectories 
will be described in this paper. 

Vehicle/Mission Description 

The X-34 Testbed Technology Demonstrator rep- 
resents a key element of the RLV program. Originally 
envisioned as an operational vehicle with the capability 
to deliver relatively small payloads (1000-2000 lb.) to 
low Earth orbit, the program was restructured in early 
1996 to focus on demonstration of RLV technologies. 6 
OSC, formerly teamed with Rockwell International in 
the initial effort, 7 was selected in June of 1996 to lead 
the redirected program. An artist’s concept of the X-34 



Fig. 1 Artist's concept of the X-34 configuration 
inflight . 
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Fig. 2 Typical X-34 flight profile. 

flight vehicle is shown in Fig. 1 . Compared to the Shuttle 
orbiter, the vehicle is relatively small with a length of 
-58 ft, a wing span of -28 ft and a height (measured 
from the bottom of the fuselage to the tip of the vertical 
tail) of - 1 2 ft. An illustration of a typical flight profile is 
shown in Fig. 2. The vehicle is designed to be air- 
launched from an L-101 1 at a Mach number of approxi- 
mately 0.7, ignite its 60 klb thrust bi-propellent (LOX/ 
kerosene) engine, accelerate to speeds up to Mach 8, 
reach altitudes up to 250 kft, and finally to coast to a 
horizontal landing on a conventional runway. 

Analysis Tools 


MINIVER 

MINIVER 8 is a versatile engineering code that uses 
various well-known approximate heating methods, to- 
gether with simplified flowfields and geometric shapes 
to model the vehicle. Post-shock and local flow proper- 
ties based on normal-shock or sharp-cone entropy con- 
ditions are determined in MINIVER through user selec- 
tion of the various shock shape and pressure options. 
The calculations can be based on perfect-gas or equilib- 
rium-air chemistry. Angle-of-attack (AOA or a) effects 
are simulated either through the use of an equivalent tan- 
gent-cone or an approximate crossflow option. 9 The flow 
can be calculated for either two- or three-dimensional 
surfaces. However, the three-dimensional effects are 
available only through the use of the Mangier transfor- 
mation 10 for flat-plate to sharp-cone conditions. MINI- 
VER has been used extensively as a preliminary design 
tool and has demonstrated excellent agreement with more 
detailed solutions for stagnation and windward acreage 
areas on a wide variety of vehicle configurations. 11 14 

LATCH 

The Langley Approximate Three-dimensional Con- 
vective Heating (LATCH) code 15 employs an inviscid/ 
boundary layer methodology that can be used to com- 
pute surface heating rates, both laminar and turbulent, 
on three-dimensional vehicles at AOA. An axisymmetric 
analog concept for three-dimensional boundary layers 
is used in conjunction with a generalized body-fitted 
coordinate system and an approximate heating method 
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developed by Zoby. 11 16,17 The LATCH code, using the 
method in Ref. 16, requires far less computational re- 
sources than a full boundary-layer solution. The LATCH 
method has previously been shown to provide accurate 
results for both wind-tunnel and flight conditions. 15,18,19,20 
An inviscid flow-field solution is used in the computa- 
tion of boundary-layer edge properties and also to pro- 
vide the surface velocities for computing the streamline 
directions. This method is described in detail in Ref. 15. 

Inviscid Solutions 

For the purpose of this study, the inviscid solutions 
required for use with the LATCH code are computed 
using either the Data-Parallel Lower-Upper Relaxation 
(DEPLUR) method 21,22 or the Langley Aerothermody- 
namic Upwind Relaxation Algorithm 23,24 (LAURA) run 
in the inviscid mode. These two codes, which take 
advantage of different computer architectures, allowed 
solutions to be run concurrently on different machines, 
generating a greater number of solutions than would have 
been possible using a single computer resource. Both 
codes utilize rectangularly- ordered structured grids as 
required for compatibility with LATCH. LATCH results 
for the peak heating condition computed using the 
DEPLUR and the LAURA inviscid solution are shown 
to be in good agreement, 19 both with each other and with 
a detailed viscous thin-layer N-S solution computed using 
LAURA. The LAURA viscous solution is described in 
detail in Ref. 25. Previous studies on a wide variety of 
vehicles including the Space Shuttle orbiter 26 , the 
Lockheed X-33 20 , HL-20 27 , and Reentry-F 8 have dem- 
onstrated the accuracy of the LAURA solutions com- 
puted for both flight and wind-tunnel conditions. 

Wind-Tunnel Test Techniques 

Advances in wind-tunnel test techniques have re- 
cently enabled the rapid acquisition of quantitative glo- 
bal heat-transfer measurements. The method, known as 
two-color relative-intensity phosphor thermography, is 
described in Refs. 5, 29, 30, and 31. Ceramic models, 
typically on the order of 12 in., can be fabricated rap- 
idly, and a wide range of parametrics, including AOA, 
sideslip, Reynolds number (ReJ, Mach number (MJ, 
and control surface deflections can be investigated. For 
the X-34 heating studies, 32 results have been obtained 
using this technique in the Langley’s 31 -Inch Mach- 10 
and 20-Inch Mach-6 air facilities. These facilities are 
described in detail by Micol in Ref. 33. Digital optical 
measurements allow the rapid processing of data using 
Merski’s IHEAT 5 code. In addition to the global images 
available, spanwise and chordwise cuts can be presented 
for comparison with engineering or detailed solutions. 

Flow visualization techniques are also used to 
complement the surface heating tests. 32 The oil-flow tech- 


nique is used to highlight surface streamline patterns. A 
stainless steel model, painted black and coated with white 
pigmented oils, was tested in both the Mach-6 and the 
Mach- 10 facilities. The Mach-6 facility is equipped with 
a schlieren system which uses the density variation in 
the flowfield to highlight the existing shock patterns. 
These techniques, although qualitative in nature, repre- 
sent an important contribution to the heat-transfer analy- 
ses by highlighting regions of the vehicle that may re- 
quire detailed analysis. 

TPS Considerations 

OSC has chosen a TPS which utilizes thermal blan- 
kets over the acreage of the vehicle and ceramic tiles in 
the stagnation regions of the nose and the leading edges. 
This decision allows OSC to take advantage of a tech- 
nology which has been demonstrated on the Space 
Shuttle orbiter, together with developments that have 
advanced these systems to improve the operational char- 
acteristics. A chief advantage of the flexible blankets is 
that they can be adhesively bonded to a non- smooth sur- 
face. Also, they should be fairly tolerant of the over tem- 
peratures that may be associated with heating uncertain- 
ties. However, the quilted-blanket design presents a fairly 
rough surface, and is more likely to induce turbulent flow 
than would a smooth surface. Also, use of blankets in 
areas of high shear stress has not been demonstrated. 
The TPS layout is shown in Fig. 3. Three types of blan- 
kets are utilized on the X-34 vehicle, depending on the 
peak temperatures expected. The First, Advanced Flex- 
ible Reusable Surface Insulation (AFRSI), is used for 
regions expected to experience temperatures no greater 
than 1500°F. This blanket is flight certified and currently 
in use on the leeward regions of the Shuttle orbiter. As 
can be seen in Fig. 3, the second type of blanket is a 
High-Temperature AFRSI (HT-AFRSI) designed for 
reuse at temperatures up to 2000°F. This blanket is used 
over the majority of the windward surface. In this blan- 

AFRSI 



Fig. 3 X-34 TPS blanket layout . 
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ket, the Astroquartz fabric and Q-felt insulation of the 
standard AFRSI are replaced with Nextel 440 and Saffil 
insulation, respectively, increasing its temperature ca- 
pability. All blankets are assumed to be coated with Pro- 
tective Ceramic Coating (PCC), a high-emissivity coat- 
ing designed to improve radiant characteristics at high 
temperatures. Nomex Felt blankets, also known as Flex- 
ible Reusable Surface Insulation (FRSI), are to be used 
in the lowest temperature regions, up to 700°F (prima- 
rily the side fuselage and the upper surface). The blan- 
kets and coating are described in the Ames TPS materi- 
als portion of the TPSX Database. 34 SIRCA 2 tiles are 
used in the high-temperature stagnation areas. Design 
and analysis for the tile regions was performed at ARC. 2 4 
The LaRC effort concentrated on the blanket areas of 
the vehicle. Areas of concern immediately noted were: 
roughness-induced turbulent heating, the design of the 
interface between the tile and blanket systems, potential 
areas of shock impingement and shock interactions, use 
of blankets in the high-dynamic pressure regions of de- 
flected control surfaces, the zone of transition of the fu- 
selage cross-section from circular to square, and the vor- 
tices which will likely emanate from the strake fuselage 
juncture. 


Overview of Aerothermal Approach 
for TPS Design 


This paper details the aerothermal environments 
prediction process as it has been applied to the X-34 ef- 
fort. In the year and a half since LaRC began this effort, 
the thermal environment predictions have evolved from 
preliminary heating estimates limited to the stagnation 
areas and the windward centerline, to environments for 
eighty representative points on the vehicle, including off- 
centerline windward, leeward, and side fuselage loca- 
tions. Prior to the acquisition of any wind-tunnel or com- 
putational data, a conservative engineering approach was 
employed. As wind-tunnel data became available, 32 com- 
parisons were made with engineering predictions and 
the engineering model was refined to reduce the conser- 
vatism in the predictions. An approach that coupled the 
experimentally-measured heating distributions over the 
body with the time-dependent engineering predictions 
along the centerline was developed. This methodology 
was used to make the first predictions of the thermal 
environments over the acreage of the vehicle. The pro- 
gram schedule and computational resources limited the 
number of detailed Navier- Stokes (N-S) heating solu- 
tions that could be used to benchmark the engineering 
methods. Instead, a coupled inviscid/boundary layer 
method 19 was used. More points along the trajectory 
could then be computed because the inviscid/boundary 
layer method requires fewer computational resources 


than a N-S method. Inviscid solutions were coupled with 
an approximate heating method 16 using the LATCH 15 
code to provide heating distributions over the surface of 
the vehicle. The inviscid solutions were chosen to bound 
the heating-rate, AOA, freestream Mach-number, and 
Reynolds-number ranges over the flight profile, includ- 
ing the peak heating condition. Finally, the engineering 
code, MINIVER, 8 was used to interpolate the LATCH- 
generated spatial distributions in Mach number, AOA, 
and Reynolds number over the time of the trajectory. 
Thus, the time- dependent thermal environments neces- 
sary for design of the flexible TPS to be used over the 
acreage of the vehicle could be predicted through an in- 
tricate integration of wind-tunnel, detailed computa- 
tional, and engineering methods. The thermal environ- 
ments generated in this manner have been used to define 
the TPS requirements for the X-34 acreage areas 3 pro- 
tected by the flexible systems. 

Preliminary Aerothermal Predictions 

The first flight of the X-34 Technology Testbed 
Demonstrator is currently scheduled to occur prior to 
the end of 1998, approximately two and one-half years 
from the date of the award of the final contract In this 
remarkably short time frame, a vehicle system must be 
designed, materials ordered, components manufactured 
(including the engine which has yet to be demonstrated), 
ground facilities built, operational issues addressed, the 
vehicle assembled, the carrier aircraft modified, the ve- 
hicle mated, and ready to fly. The extraordinary pace of 
this program required that thermal environments be gen- 
erated in parallel with design of the vehicle and devel- 
opment of the flight profile. The lead time required to 
set up for manufacture of the thermal blankets for the 
acreage and the tiles for the stagnation areas required 
preliminary definition of the TPS requirements very early 
on in the program. In light of these considerations, it 
was necessary to proceed with engineering estimates of 
the heating environments before any wind-tunnel data 
or detailed computations were available, and before a 
legitimate flight profile had been determined. A conser- 
vative engineering approach was therefore used for ini- 
tial estimates. It was assumed that the conservatism in 
the design could be reduced as the fidelity of the models 
improved with introduction of additional wind-tunnel and 
computational results. 

The MINIVER code has been shown to produce 
fairly accurate predictions of the time-dependent ther- 
mal environments for a variety of vehicle concepts, 11 14 
particularly for the windward surface regions and stag- 
nation areas. Of course, the engineering techniques used 
will not predict localized flow phenomenon such as heat- 
ing due to shock interactions or imbedded vortices. Ar- 
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eas of potential concern were identified and recom- 
mended for detailed computational study or wind-tun- 
nel testing as appropriate. The inherent roughness of the 
quilted blanket system coupled with the low-altitude, 
high-Reynolds-number flow regime led OSC to make 
the decision that the acreage TPS be designed for fully 
turbulent flow and heating levels. Langley concurred 
with this decision and, although some laminar calcula- 
tions were performed to validate computational efforts, 
the environments provided for blanket design were all 
fully turbulent. 

When simple engineering techniques are used, pre- 
diction of accurate thermal environments for a vehicle 
such as the X-34 is very much dependent on the experi- 
ence of the thermal analyst and recognition of the limi- 
tations of the codes. 11 The windward centerline 
calculations were performed assuming sharp-cone en- 
tropy conditions, flat-plate heating with the Mangier 
transformation used to correct for a three-dimensional 
boundary layer effects. When seen in cross-section, the 
windward surface of the X-34 vehicle is extremely flat. 
This, coupled with the low angle-of-attack trajectory 
profiles suggested limited crossflow effects and no 
crossflow correction was included in the initial heating 
estimates. As stated previously, all acreage environments 
were based on fully turbulent flow conditions. 

The initial trajectory used was denoted “WAG1 8D”, 
and is presented in terms of Mach number, altitude, and 
AOA in Fig. 4. A more detailed discussion of the trajec- 
tory development and analysis will be discussed in a later 
section. It should be noted however that, of all the tra- 
jectories examined, the WAG18D was found to be the 
most severe in terms of both heating rates and loads. A 
summary of the early engineering predictions is presented 
in Fig. 5. Figure 5(a) shows the radiation equilibrium 
temperature time histories for the nosecap and the lead- 
ing edges of the strake and wing. As noted previously, 
the nosecap and wing leading edge are constructed of 
SIRCA tile. An emissivity of 0.86 is therefore assumed 
for these areas. Comparison with similar calculations 35 
performed by ARC show reasonable agreement. Radia- 



Time, sec 


tion equilibrium temperatures shown throughout this 
paper are presented only as a general indicator of the 
thermal environment and do not account for the thermal 
mass of the system. Radiation equilibrium temperatures 
represent the maximum expected surface temperature 
values. Figure 5(b) shows the peak radiation equilibrium 
temperature and the total-heat-load distribution predicted 
for the windward centerline. A constant emissivity of 0.8 
is assumed. The emissivity of the PCC coating used on 
the blankets varies from about 0.75 to 0.85 over this tem- 
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a) Stagnation area temperature-time histories . 
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Fig. 4 WAG18D trajectory profile. 


Fig. 5 Preliminary aerothermal predictions . 


5 



perature range. It can be seen that, even for this severe 
trajectory case, the centerline heating environments 
should be well within the capability of the HT-AFRSI. 
However, note that heating levels in the region where 
the flow expands around the chine of the vehicle would 
be expected to exceed the centerline values considerably. 

The actual environments provided to the blanket 
designers, Oceaneering Space Systems (OSS), were pre- 
sented in terms of the heat transfer coefficient and re- 
covery enthalpy time histories for several points along 
the centerline. OSS coupled these environments with a 
thermal model based on the TPS material properties and 
design to generate the actual thermal loads and determine 
the required the blanket thicknesses. 3 Figure 5(c) illus- 
trates typical spanwise heating distributions at several x/L 
stations (where x is measured longitudinally from the tip 
of the nose, and L is the vehicle reference length) on the 
windward surface. These data were also provided to aid 
OSS in blanket design . The laminar distributions, based 
on the 0.0153-scale model 32 of the X000912 (a pre- outer- 
mold-line freeze configuration), were available for use 
in the initial blanket sizing and give some indication of 
the amplification of the heating from the centerline to the 
vehicle chine. Experience indicates that, although the 
turbulent heating levels will be higher, the amplification 
from the centerline to the chine will be reduced in turbu- 
lent flow, 36 partially due to the reduced effect of 
crossflow on turbulent heating levels. The 20-deg AOA 
case was selected because it generally showed the high- 
est amplification from the centerline to the chine. OSS 
was able to couple the centerline time histories with the 
wind-tunnel amplification to make an initial determina- 
tion of the blanket requirements over the windward sur- 
face. The initial blanket sizing and layout by OSS 3 was 
based entirely on these initial heating estimates and pre- 
liminary wind-tunnel heating distributions because no 
detailed computational solutions would be available un- 
til the outer mold line was frozen and surface and vol- 
ume grids could be generated. 

As additional wind-tunnel data became available it 
was used to confirm or suggest modification of the 
present engineering methodology. A sample of these data 
from Ref. 32 is presented in Fig. 6. Oil-flow views, simi- 
lar to that presented in Fig. 6(a), showed evidence of 
inflow at AOA of 20 deg and lower. This result sug- 
gested the use of a modification of Baranowski’s 
crossflow method 9 to account for the thickening of the 
boundary layer along the centerline due to inflow. This 
procedure has been shown to provide fairly good agree- 
ment with more detailed methods on vehicles such as 
the X-33 20 which also experiences inflow at the lower 
AOA. Schlieren views, such as that illustrated for a 15- 
deg AOA case in Fig 6(b), show a relatively strong coni- 
cal shock, thus confirming the choice of sharp-cone en- 


tropy for the engineering analysis. Figure 6(c) shows a 
typical windward centerline experimental heating dis- 
tribution normalized by a reference Fay Riddell 37 stag- 
nation point value. These data show clear evidence of 
transition at about the 0.25 x/L station. The results of 
the MINIVER heating predictions are shown to be in 
good agreement with the laminar and turbulent experi- 
mental data. The thermal environments were updated to 
reflect the minor modification in methodology. The dif- 
ferences in maximum temperatures and total heat loads 



a) Oil-flow pattern (windside view, a = 75°). 



b) Schlieren image (side view, a = 15°). 



c) Comparison of MINIVER heat-transfer results with 
experimental data (a- 20°). 

Fig. 6 Preliminary wind-tunnel results. 
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were so insignificant that no modification of the blanket 
sizing was required. Although, wind-tunnel tests often 
cannot match the actual flight Reynolds number condi- 
tions or flow chemistry, comparisons such as that shown 
in Fig. 6(c) validate the analysis tools for application to 
actual flight conditions. 

Trajectory Considerations 

The rapid pace of the X-34 program required the 
environment predictions to be updated repeatedly as the 
vehicle configuration and trajectory were developed. All 
trajectories discussed herein have been generated by 
OSC. The initial trajectory (WAG1 8D), available in July 
of 1996 and shown in Fig. 4, was based on a 3-degree- 
of-freedom analysis (control deflections and trim require- 
ments not necessarily considered) using the vehicle 
weight and performance available at that time. The maxi- 
mum Mach number is about 9.3 and the maximum AO A 
is 25 deg. An updated trajectory profile, designated the 
DRM2, was made available at the Systems Requirements 
Review early in the fall of 1996. This trajectory, illus- 
trated in Fig 7(a), reached a Mach number of only about 
7.4 with the same maximum AO A (25 deg), apparently 
indicating either reduced performance of the vehicle or 
a change in mission requirements. One feature which 
should be noted in each of these early trajectories is the 
“bounce” that occurs at about 1 10 to 120 kft. Once the 
outer mold line was “frozen” in December of 1996, 
higher fidelity trajectories could be developed. As part 
of this effort, a guidance system would need to be de- 
signed. OSC elected to eliminate the “bounce” in the 
trajectory in order to reduce the demands on the guid- 
ance system. The two trajectories shown in Figs. 7(b) 
and 7(c) and are designated the XI 004601 and the 
XI 004701 trajectories respectively. The altitude-time 
histories clearly show that the “bounce” has been suc- 
cessfully eliminated from both. Similar, to the DRM2, 
the XI 004601 has a peak Mach number of 7.4 and a 
maximum AO A of 25 deg. The XI 004701, on the other 
hand, has peak Mach number of 8.6 and maximum AOA 
of 30 deg. It would be expected that the latter trajectory 
would yield the more severe windward heating environ- 
ment with a possible reduction on the leeward side due 
to the increased AOA. However, at the request of OSC, 
LaRC provided all thermal environments based upon the 
X1004601 flight conditions. Possible reasons for OSC’s 
initial selection of this trajectory to be used for blanket 
design will be discussed later in this section. 

Figure 8 illustrates a comparison of the four trajec- 
tories in terms of reference heating levels. Reference 
values are based on a Fay-Riddell 37 calculation of the 
stagnation heating to a 1 -ft sphere, corrected for a hot 
wall. These results indicate that the original trajectory 


(WAG18D) would probably subject the vehicle to the 
highest temperatures, requiring more extensive use of 
the high temperature blankets, while the XI 004701 
would result in the greatest heat loads, possibly requir- 
ing greater TPS thicknesses. It should also be noted that 
the peak heating levels for all the newer trajectories are 
similar, and significantly lower than those for the initial 
design trajectory. Hot wall heating results presented in 
Fig. 9 show that the XI 004701 trajectory presents some- 
what more severe environments than the XI 004601. A 



a) DRM2. 



b) X1004601 (TPS blanket design trajectory). 



c) X1004701. 

Fig. 7 Trajectory profiles. 
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Fig. 8 Reference heating-rate comparisons for the 
design trajectories. 

cold wall heating comparison, presented in Fig. 9(b) for 
the same trajectories, would suggest the opposite con- 
clusion. Trajectory codes such as the Program to Opti- 
mize Simulated Trajectories (POST) 38 typically use a 
cold-wall heating indicator and it is likely that such cold 



wall results led to OSC’s decision to design the blankets 
based on the XI 004601 trajectory. For flight conditions 
similar to those for the X-34 where the ratio of the wall 
enthalpy to recovery enthalpy is fairly high (~0.4), care 
must be exercised if cold wall heating rates are to be 
used as a measure of heating levels during trajectory 
development. 

Although the acreage areas are to be designed for 
fully turbulent flow, the trends seen in the reference heat- 
ing levels will hold in these regions as well This fact is 
evident in Figs. 9(c) and 9(d) which illustrate the peak 
temperature and heat-load distributions respectively 
along the windward centerline. The temperatures and the 
loads are generally greater for the XI 004701 trajectory 
than for the X1004601. 

Case Selection for Benchmark Solution 

All flight conditions selected for detailed analysis 
were chosen along the XI 004601 trajectory profile at 
the request of OSC. Computational resources limited the 
number trajectory points to a total of six. The trajectory 



a) Hot-wall reference heating-rate time histories. c) Windward-centerline peak-temperature 

distributions. 




b) Cold-wall reference heating-rate time histories. d) Windward-centerline peak-heat-load distributions. 


Fig. 9 Heating comparison of the XI 004601 and the X 1004701 trajectories. 
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points, the case numbers, and the associated parameters 
are listed in Table 1 . Note, Case 2 was also calculated at 
the peak heating condition (t = 330 sec), but is not shown 
as it was not used to anchor the engineering results. 
LATCH solutions were available at all conditions, 

Table 1 LATCH solutions 


whereas viscous solutions were only available for Cases 
1 and 5. Case 8 will be discussed later in the paper and 
was not used in the initial analysis. Cases were chosen 
to include the peak heating on ascent and descent, as 
well critical AOA, and an appropriate range of Reynolds 
numbers. Figure 10 shows the critical trajectory param- 
eters as functions of time. Included are; altitude- 10(a), 


Time 

(sec) 

Case # 

AOA 

(deg) 

Mach 

number 

Re/ft 

q ref ( hot ) 

(Btu/ft 2 -sec) 

145 

7 

9 

6 

1.9 x 10 4 

3.0 

152 

4 

11 

6.83 

1.4 x 10 4 

3.1 

330 

1* 

23 

6.32 

2.8 x 10 5 

8.3 

340 

5* 

15.22 

6 

3.6 x 10 5 

7.3 

355 

3 

8 

5.8 

3.8 x 10 5 

6.5 

578 

6 

6.48 

3.6 

7.5 x 10 5 

0.9 

334 

00 

* 

* 

20.3 

6.17 

3.2 x 10 5 

7.8 


* Viscous Solutions Available 
** Additional LATCH solution obtained to check 
interpolation validity 



a) Altitude. 



b) Mach number. 



c) ol 



d) Re^/ft. 



e) Hot-wall reference heating rate . 


Fig. 10 X1004601 trajectory profile. 


9 



Mach number- 10(b), AOA- 10(c), Reynolds number/ft- 
10(d), and reference heating rate- 10(e). The LATCH 
(inviscid/boundary layer) solutions are identified on each 
plot, as are the points at which the viscous solutions were 
generated. Although the solutions appear to be 
“clumped” near peak heating, examination of the AOA 
profile indicates the reason for this choice. Accurate cir- 
cumferential heat transfer distributions around the body 
are critical to the success of the methodology employed 
in this paper. Because the majority of the high-heating 
portion of the trajectory does not vary widely in Mach 
number and the flow state is defined to be turbulent, the 
chief influence on the heat transfer distributions is the 
AOA. As shown in Fig. 10(c), the AOA range is cap- 
tured fairly well over the high heating portion of the entry 
(see Fig. 10(e)). The point at 330 sec (23 deg) should 
adequately represent the peak heating to the windward 
surface while that at 355 sec (8 deg) should represent 
the highest heating on the leeward surface due to the 
low AOA. The case at 340 sec (15 deg) provides heat- 
ing distribution definition where the AOA is changing 
rapidly. It is noted that the actual peak heating condition 
is not calculated as the result of the apparent selection 
of the peak heating point based on the cold wall values. 
Comparing Figs. 9(a) and 9(b) it is apparent that the peak 
in hot wall heating occurs prior to that for the cold wall, 
and therefore the points chosen fall on the downward 
slope of the reference heating profile. The lack of a so- 
lution exactly at peak heating should not impact the time- 
dependent heating significantly, as the engineering so- 
lutions are run at the actual flight conditions with the 
heating distributions about the body based on the LATCH 
solutions. The difference in distributions at 25 deg AOA 
compared to 23 deg is expected to be small. 

Detailed Aerothermal Predictions 
Code Validation 

Comparisons are presented in this section to dem- 
onstrate the validity of the more detailed techniques, as 
well as to further demonstrate the applicability of the 
MINIVER code. For the comparisons with the ground- 
test experimental heating data, several figures are used 
for the sake of clarity. Initially, the data are compared to 
the MINIVER predictions. Then the MINIVER results 
are compared to the LAURA and to the LATCH results 
respectively. Finally, comparisons of the predicted re- 
sults are made at a flight-test condition. 

The aerothermal environments used for the initial 
design of the blanket TPS were based strictly on engi- 
neering estimates coupled with preliminary wind-tunnel 
distributions. As the study progressed and the vehicle 
geometry was defined, wind-tunnel test data, LATCH 
solutions, and viscous solutions became available for 



a ) Experimental data ( thermographic phosphor). 



b) LAURA viscous results . 



c) LATCH results 

Fig. 11 MINIVER results compared at wind-tunnel 
test conditions. 
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comparison to the engineering predictions. All 
engineering predictions shown here are made using the 
updated methodology as described in the Preliminary 
Aerothermal Predictions section. In Fig. 11 the 
MINIVER engineering predictions for the windward 
centerline of a 0.183 scale model of the X-34 config- 
uration X0001215 at a wind-tunnel condition of Mach 
6, an AOA of 15 deg, and unit Reynolds number of 7.9 
million are compared to thermographic phosphor data 32 
and analytical predictions. 19,25 In Fig. 1 1(a) the MINIVER 
predictions are compared with the thermographic phos- 
phor data where transition is evident at an x/L of -0.25. 
Both laminar and turbulent MINIVER solutions are 
shown. Once again the data are presented in the form of 
h/h^ (where h is the local heat transfer coefficient and 
h F is the Fay-Riddell reference value). Good agreement 
with the data is evident for both the laminar and turbulent 
results. MINIVER’s overprediction of the laminar results 
in the forward portion of the body is expected and results 
primarily from the use of sharp-cone entropy conditions 
in this region where the streamlines are more likely to 
have passed through a higher angle shock near the nose 
of the vehicle. In Fig. 1 1(b) the MINIVER calculations 
are compared with the LAURA solutions. Once again 
good agreement is noted, particularly downstream where 
sharp cone entropy is most appropriate. Finally, in Fig 
1 1(c), the MINIVER predictions are compared with those 
from the inviscid-boundary layer solution of LATCH. 
Reasonable agreement is demonstrated once again, 
although it is apparent by comparison with Fig. 11(a) 
that the LATCH solutions tend to over-predict the data 
and the viscous solutions at centerline stations down- 
stream of approximately 0.6 x/L. The difference appears, 
however, to be less than 1 0% for the turbulent data which 
is of primary interest here. 

Figure 12 presents the MINIVER engineering pre- 
diction of turbulent heating rates along the windward 
centerline for the full scale vehicle at the peak-heating 
flight condition (t = 330 sec) on the XI 004601 trajec- 
tory. Downstream of the 200 in. station where the sharp- 
cone-entropy assumption is more appropriate, good 
agreement is shown with the viscous LAURA solution 
and the LATCH results computed for the same condi- 
tion. Similar levels of agreement are shown for the ra- 
diation equilibrium temperatures not shown here. 

A MINIVER generated temperature-time history is 
shown in Fig. 13 for a representative location on the 
windward centerline for the XI 004601 trajectory (used 
for blanket design). The LATCH-generated values are 
shown for comparison. It is evident, that, even without 
adjustment, the MINIVER predictions of the tempera- 
tures is in excellent agreement with the LATCH predic- 
tions. Where any noticeable difference exists MINIVER 
tends to overpredict, resulting in the conservative esti- 



Fig. 12 MINIVER/LAURA/LATCH predictions 
compared at flight condition ( Case 1 ). 



Time, sec 

Fig. 13 MINIVER temperature comparison with 
LATCH results . 

mate desired. It should be noted that, particularly in re- 
gions of steep temperature gradients, a certain degree of 
error is introduced because the LATCH solutions are 
allowed to come to equilibrium at a constant wall tem- 
perature whereas the MINIVER wall temperature lags 
the radiation equilibrium value slightly. This case is typi- 
cal of the agreement shown along the windward 
centerline for stations downstream of the nose region. 

Comparisons of the MINIVER predictions with 
wind-tunnel data and detailed solutions at flight condi- 
tions suggest that the initial heating predictions provided 
for blanket design will be somewhat conservative over 
the forward portion of the vehicle, but predict the envi- 
ronments over the remainder of the acreage with rea- 
sonable accuracy. This, in combination with the more 
severe trajectory originally considered, indicates that the 
blanket TPS designed using the original environments 
provided and the WAG18D trajectory should be suffi- 
cient to handle both the maximum temperatures and to- 
tal heat loads associated with the refined trajectories now 
available. Demonstration of the validity of the engineer- 
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ing predictions is critical to early definition of accurate 
thermal environments. However, once a sufficient num- 
ber of detailed solutions become available, the accuracy 
of the engineering predictions is less important as they 
are anchored to the detailed solutions and the engineer- 
ing methods are then used primarily as an interpolation 
tool for the trajectory conditions which fall between those 
for which detailed solutions are available. 

Methodology 

Limited computational resources together with the 
time constraints imposed by the program schedule pre- 
cluded the generation either of sufficient number of 
LAURA viscous solutions or a sufficient number of 
LATCH solutions to provide the time-dependent ther- 
mal environments necessary for TPS design. Instead, two 
LAURA solutions 25 were generated at or near peak heat- 
ing for the appropriate AOA and flight conditions. 
LATCH solutions were computed for the same cases and 
at four additional flight conditions. The LATCH solu- 
tions are shown to compare well with the corresponding 
viscous cases in Ref. 19. The six LATCH cases were 
then used to anchor the engineering solutions. Finally, 
the MINIVER solutions, thus anchored to the LATCH 
solutions, were used to interpolate in time, Mach, 
Reynolds number, and AOA to provide the required heat- 
ing time histories at eighty points on the body. A de- 
scription of the way in which this process was imple- 
mented follows. 

Centerline Adjustments: Figure 14 shows a typical 
comparison of MINIVER and LATCH heat-transfer co- 
efficients for the windward centerline. This particular 
figure shows the heating comparison at 330 sec on the 
XI 004601 trajectory. Plots such as this are used to de- 
termine the factors necessary to correct the MINIVER 
prediction to match the LATCH value. The ratio of h^ 
(the LATCH value) to h~" (the MINIVER prediction) is 
determined for approximately twenty locations along the 
windward centerline from a point immediately aft of the 
nosecap, -15 in., to the end of the LATCH solutions at 
approximately 560 in., including the point identified at 
an axial station of 350 in. This procedure is repeated for 
each condition for which a LATCH solution is avail- 
able, a total of six in this case. 

Off-Centerline Adjustments: Once the adjustment fac- 
tors have been determined for the centerline, an addi- 
tional correction must be determined to adjust the pre- 
dictions to match the heating distributions at locations 
away from the centerline. Figure 15 presents typical cir- 
cumferential heating distributions for three different 
times along the XI 004601 design trajectory represent- 
ing AOA of 8.0, 15.2, and 23 deg. The data are pre- 



Fig. 14 MINIVER/LATCH windward-centerline heat- 
transfer coefficient comparison for peak heating case 
(Case 1). 



Fig. 15 Typical turbulent circumferential heating 
distributions from LATCH. 

sented as an h/h CL distribution for a cut plane 19 located at 
an axial station of 350 in., corresponding to the centerline 
location identified in the previous figure . The increas- 
ing influence of inflow with reduced AOA is evidenced 
by the corresponding increase in heating amplification 
at the most outboard locations, (see Ref. 19 for other 
sample heating distributions and comparisons with the 
viscous results) For each solution and each cut plane, 
factors were evaluated at a sufficient number of spanwise 
locations to adequately define the distribution (typically 
three to five points) such that linear interpolation be- 
tween points would yield reasonable agreement with the 
detailed solution. A point representing the greatest am- 
plification factor at the peak heating condition, 330 sec, 
was always chosen as one of the points in order to en- 
sure that the peak temperature at or near the leading edge 
would be captured for the blanket TPS design. Once the 
h/h CL factors have been determined for each point and 
each solution, these were multiplied by the original 
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centerline factor (h CL /h min ) to determine the required 
MINIVER adjustment factor. The local heat transfer co- 
efficient (h) is calculated using the following equation: 

h = h x(h„_/h ) x (h/h ) (1) 

min v CL min ' x CL 7 v 7 

where: h^/h . is evaluated at the appropriate 
centerline location 

h/h cl is evaluated at the desired spanwise 
position along the cut plane 
h is calculated as a function of time using 
standard engineering techniques. 

Thus a MINTVER-generated centerline heating-rate time- 
history can be corrected to yield the time-dependent ther- 
mal environment for any point on the body for which 
the appropriate set of factors has been determined. The 
next session discusses the appropriate application of these 
factors given a limited number of detailed solutions so 
that a valid interpolation of the heat transfer can be per- 
formed. 

Solution Application: Ideally a database which con- 
tains solutions for a range of AO A, Mach number, and 
Reynolds number would be created. However, available 
time and resources precluded that approach here and 
limited the data set to six solutions at conditions on the 
X1004601 trajectory (see Table 1). Engineering judg- 
ment was required to choose the appropriate application 
of this limited solution set over the trajectory. TPS sizing 
required computations be performed at small time incre- 
ments over the entire time of the trajectory (-700 seconds 
from release of the vehicle through ascent to 264 kft and 
descent back to Mach -1.5). Referring to Fig. 10(e), it is 
apparent that the majority of the heating will occur be- 
tween 100 and 600 sec making the case selection before 
and after those times less critical. Wherever a solution 
was available it was applied. Table 2 lists the trajectory 
times at which the six solution cases were applied. Addi- 
tionally, this table lists the corresponding range of AOA, 



Fig. 16 Typical temperature -time history for centerline 
adjusted to windward strake location. 

Mach, and Reynolds numbers. In the time from 100 to 
600 sec, AOA was given priority in case selection. Com- 
parison of the conditions in Table 2 with those in Table 
1 for the design cases shows that Mach numbers gener- 
ally did not vary significantly from the design point used. 
In some cases Reynolds number is seen to vary substan- 
tially from the design case. However, because a fully 
turbulent boundary layer was assumed for all heating 
distributions this variation was not considered critical. 
Where two times are listed, the distributions were held 
constant (but not the MINIVER reference). Times were 
generally chosen such that the flight AOA corresponded 
to those of the detailed solutions, thus allowing linear 
interpolation from one time step to the next. Figure 16 
shows a typical radiation equilibrium temperature pro- 
file resulting from this process. The dashed line repre- 
sents the MINIVER centerline time history without ad- 
justment, while the solid line represents the heating asso- 
ciated with an outboard point at the corresponding cut 
plane. This particular location represents a windward 
point near the edge of the strake. The points on the body 
at which detailed time histories were provided for blan- 
ket design will be discussed in a later section. 


Table 2 Trajectory application of LATCH solutions for interpolation 


Time 

(sec) 

Case # 

AOA 

(deg) 

Mach 

number 

Re/ft q ref range (hot) 

(Btu/ft 2 -sec) 

0-137 

6 

4-14 

.7-5.2 

3.0 x 10M.5 x 10 6 

.1-2.5 

145 

7 

9 

6 

1.9 x 10 4 

3.0 

150-240 

4 

10.5-11.1 

6.7-7. 3 

1.1 x 10M.5 x 10 4 

.9-3.2 

242-331 

1 

23-25 

6.3-7. 3 

1.1 x 10 3 -2.8 x 10 5 

.9-8.3 

340 

5 

15.22 

6 

3.6 x 10 5 

7.3 

355 

3 

8 

5.8 

3.8 x 10 5 

6.5 

578-714 

6 

6. 6-5. 5 

1.6-3.5 

7.7 x 10 5 -8.5 x 10 5 

0-.9 
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Fig . 17 Typical heat-transfer-coefficient and recovery- 
enthalpy time histories for blanket design. 

Aerothermal Environments: The aerothermal envi- 
ronments for blanket design were provided in the form 
of time-dependent heat-transfer coefficients and recov- 
ery enthalpies (adiabatic wall enthalpies) similar to those 
illustrated in Fig. 17 for a typical centerline location. 
The blanket designers were then able to perform a ther- 
mal balance, coupling the applied environments provided 
with the material response of the TPS, to determine the 
blanket thickness necessary to ensure that structural tem- 
peratures limits would not be exceeded. 3 The detailed 
solutions 19,25 were used to determine the distribution of 
the different blanket materials (HT-AFRSI, AFRSI, and 
FRSI) depending upon the maximum surface tempera- 
tures anticipated. 

As noted previously, radiation equilibrium tempera- 
tures are presented here only as a general indicator of 
the thermal environment and do not account for the ther- 
mal mass of the system. In the following Results section 
of the paper, several radiation-equilibrium time histo- 
ries will be presented. The locations for the calculations 
include the wind and lee surfaces as well as the side of 
the vehicle. The results are presented to give a general 
understanding of the temperature environments for the 
X-34 blanket design trajectory (XI 004601). 

Results 

Analysis Locations 

Figure 18 shows the locations (standard body-ori- 
ented x, y, z coordinate system with the x-axis origin at 
the nose, the y axis oriented out the starboard wing, and 
z defined by the right-hand rule) selected for calculation 
of the time-dependent aerothermal environments. The 
locations are superimposed on windward, leeward, and 
side views of the X-34 vehicle. LATCH results 19 are in- 
dicated in terms of temperature contours. Fifty windside 
locations were chosen based on the peak-heating case. 
Case 1 (Mach = 6.32, AOA = 23 deg, time = 330 sec) 
and are shown in Fig. 18(a). They include locations on 


= 6.32, a = 23°, 8 ei = 0°, Turbulent, € = 0.8 



a) Windside view (Case 1). 


= 5.8, a = 0°, 5 ei = -10°, Turbulent, £ = 0.8 




b) Lee side view ( Case 3 ). 


= 6, a = 15.22°, 5 el = 0°, Turbulent, £ = 0.8 



x, in. 

c) Side view (Case 5). 

Fig. 18 Body point locations for time-dependent 
environment calculations (LATCH temperature 
contours (°F)) 

fourteen cut planes. Points along the cut planes were se- 
lected such that linear interpolation between the points 
would generally yield a fair representation of the results. 
Two additional points on the centerline of the undeflected 
bodyflap were calculated (x = 630 in. and x = 660 in.), 
but not shown because the LATCH solution did not ex- 
tend to that region. 
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Comparison of all the LATCH solutions confirmed 
that Case 3 (Mach = 5.8, AO A = 8 deg, time = 355 sec) 
represented the worst-case heating on the leeward sur- 
face, primarily due to the low AOA at relatively high 
reference-heating levels. Shown in Fig. 18(b), superim- 
posed on the LATCH results for Case 3, are the loca- 
tions of the twenty-two leeward points for which time- 
histories were provided. These included eight points on 
the wing, two on the strake and the remainder on the 
fuselage. No points were provided for the deflected 
elevons, apparent in the Fig. 18(b) at roughly the 530 in. 
station. 

Time-dependent environments at additional loca- 
tions on the side of the fuselage were also provided. 
These locations are marked on Fig. 18(c) which also 
shows the LATCH temperature contours for Case 5 
(Mach = 6, AOA = 15.22 deg, time = 340 sec). Of the 
seven side-fuselage locations calculated, an important 
one to note is that at the x location of roughly 400 in. 
and the z location of roughly 50 in. This represents the 


peak heating on this surface, and occurs at the juncture 
where the rounded forward fuselage transitions to the 
square cross-section. 

Windward Surface Heating 

Windside time-dependent aerothermal environments 
have been calculated for a total of fifty-two points. Typi- 
cal results, presented in the form of radiation equilib- 
rium time histories, are shown in Fig. 19. Figures 19(a-d) 
represent calculations at four representative cut planes; 
x = 15,70,350, and 450 in. As can be seen in Fig. 18(a), 
the forward two cut planes lie entirely on the fuselage, 
while the latter two include points on the strake and wing 
respectively. The trends in the time histories are gener- 
ally consistent, the outboard locations tracking the 
centerline results, but with some amplification in heat- 
ing. This amplification is particularly evident on the 450 
in cut plane near the leading edge of the wing (Fig. 1 9(d)). 
Assuming the radiation equilibrium temperatures pro- 
vide a reasonable representation of the surface heating 


Spanwlse Spanwise 

location, location, 




a) Cut plane - 15 in. 


c) Cut plane - 350 in. 


Spanwise 

location, 



Radiation 

equilibrium 

temperature, 

°F 


Spanwise 

location, 



b) Cut plane - 70 in. d) Cut plane - 450 in. 

Fig. 19 Typical windside temperature -time histories. 
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levels, the temperature capability of the HT-AFRSI se- 
lected for use over this region (see Fig. 3) does not ap- 
pear to be exceeded for any of the points calculated. 

Leeward Surface Heating 

Time-dependent aerothermal environments have 
been calculated for a total of twenty-two points on the 
leeside of the vehicle. Figure 20 illustrates typical time- 
dependent results for this surface at four cut planes; 
x = 15, 70, 350, and 450 in. Once again, the forward 
two cut planes lie entirely on the fuselage, while the lat- 
ter two include points on the strake and wing respec- 
tively. As expected, maximum temperatures are shown 
to be significantly lower on the lee surface than those 
seen previously on the windward surface (see Fig 19). 
Peak heating on the leeward surfaces appears to be de- 
layed slightly from that evident on the windward sur- 
face. Unlike the windward surface, heating on the leeside 
lags the peak reference heating slightly and occurs as 
the AOA drops rapidly and the leeside starts to see more 
flow. The character of the heating on the leeward sur- 
face also differs from that on the windward surface in 


that the heating, although not as high in magnitude, re- 
mains at higher levels for a longer period of time as the 
vehicle AOA is decreased (i.e. the heating does not drop 
off as rapidly after peak heating). As can be seen by com- 
paring the blanket layout in Fig. 3 with the analysis lo- 
cations shown in Fig. 18b, all these locations (with the 
exception of the point on the 450-in. cut plane near the 
leading edge of the wing (SIRCA) and the centerline 
locations on the 350- and the 450-in. cut planes (FRSI)) 
are protected with the standard AFRSI which has a maxi- 
mum use temperature of 1 500 °F. Examining the time 
histories in Fig. 20, it appears that this temperature is 
unlikely to be exceeded for any of these locations. Even 
the point near the leading edge of the wing on the 450- 
in. cut plane exceeds 1500 °F only slightly (see Fig. 
20(d)). The SIRCA tile at the leading edge should be 
sufficient to handle this thermal environment. 

The majority of the lee surface of the wing is to be 
protected by FRSI as is the aft portion of the leeward 
fuselage. Figure 21 illustrates the radiation equilibrium 
time histories for two points on the leeward surface of 
the wing. Both points appear to exceed the FRSI capa- 



Racfiation 

equilibrium 



a) Cut plane - 15 in. 


c) Cut plane - 350 in. 


Radiation 

equilibrium 

temperature, 

°F 




b ) Cut plane - 70 in. d) Cut plane - 450 in. 

Fig. 20 Typical leeside temperature -time histories. 


16 





Time, sec Time, sec 


Fig . 21 Typical leeward wing temperature-time 
histories. 


Fig. 22 Typical side surface temperature-time 
histories. 


bility of 700 °F. However, the blanket layout in Fig. 3 is 
not sufficiently detailed to determine whether these 
points (shown in Fig. 18(b)) would actually lie within 
the FRSI area. The centerline locations on the 350- and 
the 450-in cut plane (mentioned previously) do appear 
to lie within in the region protected by FRSI. The pre- 
dictions shown in Figs. 20(c) and 20(d) respectively for 
these points show maximum temperatures just within 
the FRSI capability. 

Side Surface Heating 

Time-dependent aerothermal environments have 
also been calculated for seven side surface locations, 
including a point at the 395-in cut plane where the body 
cross-section transitions from circular to rectangular. 
Figure 22 illustrates some typical results for the side fu- 
selage. Except for the fuselage transition zone, the peak 
temperatures along the fuselage fall in the 750-850°F 
range. Comparing that result with the blanket distribu- 
tion shown in Fig. 3, it is evident that the regions aft of 
300 in. which are protected by FRSI are likely to exceed 
its temperature capability of 700 °F. The AFRSI used at 
the transition zone would be sufficient to withstand the 
maximum temperatures (-1030 °F) shown in Fig. 22. 

Interpolated Result Verification 

An additional LATCH inviscid-boundary layer 
solution was run in order to test the validity of the interpo- 
lation scheme used to establish the time-dependent 
aerothermal environments shown here. The test case, 
identified as Case 8 in Table 1, was run for the flight 
condition at 334 sec on the XI 004601 trajectory. An 
inviscid solution at a Mach number of 6. 1 7 and an AOA 
of 20.3 degrees was used as input to the LATCH code to 
establish turbulent heating levels over the vehicle. Figure 
23 shows some representative heating time histories as 
predicted using the techniques described in this paper. 
The predictions are shown for three points along the 


350-in. cut plane, including the windward and leeward 
centerline and a point near the leading edge of the strake. 
The predicted heating rates and radiation equilibrium 
temperatures are shown in Figs. 23(a) and 23(b) 
respectively. The results of the new LATCH solution at 
334 sec are also shown. Errors in heating rates of less 
than 10 percent are noted. The corresponding errors in 
temperature are significantly lower, less than 2 percent. 
Heat-transfer-coefficient errors (not shown) were noted 
to be less than 3%. Similar results were observed for all 
other points examined. As noted previously, a portion 
of the error is introduced simply because the radiation 
equilibrium temperature computed along a trajectory will 
not be equivalent to one simply computed at a constant 
freestream condition. The interpolated results are shown 
to compare well with the LATCH solution generated at 
an additional flight condition, demonstrating the validity 
of the approach used in this study. Thus, the techniques 
employed here appear to enable prediction of time- 
dependent aerothermal environments with a sufficient 
degree of accuracy for TPS design. 


This paper has described the methodology by which 
the aerothermal environments for the X-34 have been 
predicted in sufficient detail to allow design of the TPS 
such that the survivability, as well as the reusability, of 
the flight vehicle is ensured to a high degree of certainty. 
Initial estimates of the aerothermal environments, pro- 
vided before any wind-tunnel or detailed computational 
data were available, are shown to be conservative. Engi- 
neering results, validated against both experimental data 
and viscous N-S solutions, demonstrate good agreement. 
A combination of engineering methods, coupled with 
inviscid CFD solutions is used to generate the detailed 
time-dependent environments over the surface of the 
vehicle. The interpolated engineering values are shown 
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a) Heating rates. 

miniver 

Predictions Error 

y. In. 

55.3 windward -0.4% 

E = 0.8 0(4) windward -1.8% 



b) Temperatures. 

Fig. 23 Verification of MINIVER interpolation 
process. 

to be in good agreement with the more detailed results. 
However, it is noted that the selection of the cases for 
detailed solutions is quite critical to the success of the 
interpolation. In particular, the AOA and Mach-number 
conditions must be matched adequately and the condi- 
tions of peak windward and leeward surface heating must 
be captured. Selection of the peak heating conditions 
based on a cold wall value has been demonstrated to be 
inappropriate for vehicles such as the X-34 that operate 
in flight regimes where the recovery to wall enthalpy 
ratio is significant. 

The results presented have been used to design and 
size the blanket TPS on the X-34. The method described 
was formulated to take optimal advantage of the strengths 
of the engineering and experimental techniques together 
with those of the detailed flowfield codes. Minimal com- 
putational resources and time were required to provide 
sufficiently accurate, time-dependent aerothermal envi- 


ronments to design the X-34 blanket TPS. This method 
could easily be automated and applied to many of the 
fast-paced programs which are typical in today’s design 
environment. In doing so, the conservatism inherent to 
any TPS designed using engineering methods alone could 
be reduced, potentially decreasing the weight and im- 
proving the payload capability. 
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